Excision boundary conditions for black hole initial data 
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We define and extensively test a set of boundary conditions that can be applied at black hole 
excision surfaces when the Hamiltonian and momentum constraints of general relativity are solved 
within the conformal thin-sandwich formalism. These boundary conditions have been designed to 
result in black holes that are in quasiequilibrium and are completely general in the sense that they can 
be applied with any conformal three-geometry and slicing condition. Furthermore, we show that they 
retain precisely the freedom to specify an arbitrary spin on each black hole. Interestingly, we have 
been unable to find a boundary condition on the lapse that can be derived from a quasiequilibrium 
condition. Rather, we find evidence that the lapse boundary condition is part of the initial temporal 
gauge choice. To test these boundary conditions, we have extensively explored the case of a single 
black hole and the case of a binary system of equal-mass black holes, including the computation of 
quasi-circular orbits and the determination of the inner-most stable circular orbit. Our tests show 
that the boundary conditions work well. 

PACS numbers: 04.20.-q, 04.25. Dm, 04.70.Bw, 97.80.-d 
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I. INTRODUCTION 

The simulation of black-hole systems necessarily starts 
with the specification of initial data. In order for such 
simulations to yield astrophysically relevant results, the 
initial data must be constructed to be astrophysically re- 
alistic. Achieving this is the goal of efforts being made 
to improve black-hole, and in particular black-hole bi- 
nary, initial data. It has become clear that all of the 
freely specifiable pieces of the initial data, including the 
boundary conditions, must be chosen carefully to respect 
the physical content of any system we wish to simulate. 
In this paper we will focus on the boundary conditions 
that are required when a black hole's interior is excised 
from the initial-data domain. 

In Ref . , one of the authors proposed a set of bound- 
ary conditions that were intended to yield a black hole 
that was in quasiequilibrium. These conditions were cho- 
sen to be consistent with the desire to create a binary 
system fully in quasiequilibrium. It is reasonable to ex- 
pect such a system will be astrophysically realistic if the 
black holes in the binary are sufficiently far apart and 
in a nearly circular orbit. In this paper, we refine and 
extensively test these boundary conditions in the cases of 
a single black hole and a pair of equal-mass black holes 
in a binary system. 

The most significant refinement of the quasiequilib- 
rium boundary conditions over the original version in 
Ref. is to the procedure for specifying the spin of each 
black hole. The analysis below shows that the spin must 



'Electronic address: cookgb@wfu.edu 



^Electronic address: harald@tapir.caltech.edu 



be chosen in a very specific way in order to be compatible 
with the assumptions of quasiequilibrium. Fortunately, 
the procedure still allows for a completely arbitrary spec- 
ification of the spin and this procedure is compatible with 
any choice of the conformal three-geometry. 

A significant result from our tests on the original set 
of quasiequilibrium boundary conditions is that the pro- 
posed lapse boundary condition is not viable. We will 
show below that this boundary condition is degenerate 
when combined with the other quasiequilibrium bound- 
ary conditions. Furthermore, the nature of this degen- 
eracy can be easily understood. We conjecture that the 
boundary condition on the lapse is not fixed by quasiequi- 
librium considerations but is, rather, a part of the initial 
temporal gauge choice. Below, we provide analytical and 
numerical evidence to support this conjecture. 

For a single black hole, the quasiequilibrium bound- 
ary conditions allow for the construction of initial data 
that yield true stationary spacetimes. Doing so, however, 
requires that all of the freely specifiable data be cho- 
sen in a way that is compatible with a stationary black 
hole. In particular, it requires that the conformal three- 
geometry be chosen correctly. Unfortunately, there is 
still no general prescription for choosing an appropriate 
conformal three-geometry. Because of this, we have cho- 
sen to perform all of our tests on a flat conformal three- 
geometry. By testing the boundary conditions for the 
three cases of a single static, spinning, and boosted black 
hole, we will be able to test both the quasiequilibrium 
boundary conditions and the effect that the assumption 
of conformal flatness has on the resulting initial data. 
For equal-mass black hole binaries, we will extensively 
test the special cases of corotating and irrotational black 
holes. The numerical results we obtain will be compared 
with post-Newtonian results and previous numerical re- 
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suits for both cases. 

We begin in Sec^with a review of the conformal thin- 
sandwich decomposition of the constraints, and then de- 
rive the quasiequilibrium boundary conditions in Sec. lIIII 
In Sees. II Vl and IVl we apply the boundary conditions to 
the cases of a single black hole and to equal-mass black 
hole binary systems. Finally, in Sec. IVII we further ex- 
plore the effectiveness of these boundary conditions. 



II. THE CONFORMAL THIN-SANDWICH 
DECOMPOSITION 

In this work, we will use the standard 3+1 decompo- 
sition with the interval written as 



ds 2 



-a 2 dt 2 + •y ij {dx i + {3 l dt){dx J + ftdt), (1) 



where 7y is the 3-metric induced ona( = const, spatial 
hypersurface, a is the lapse function, and /3 l is the shift 
vector. The extrinsic curvature of the spatial slice, Kij, 
is defined by 



K„ 



1 <5 P r 



(2) 



where C n denotes the Lie derivative along the unit nor- 
mal to the spatial slice, n^. Einstein's equations, in vac- 
uum, then reduce to four sets of equations. Two are 
evolution equations for the spatial metric and extrinsic 
curvature: 



= -2aKij + 2%p } 



(3) 



and 



dtlCy = -%V j a + a[R ij -2K ii K? + KK ij ] 

+ f3 l %K l3 +2K i(l Vj)(3 l - (4) 
The remaining two are the constraint equations 



R + K 2 - K i:j K ij = (5) 



and 



Vj(K ij -j ij K) =0. 



(6) 



Here, Vi, Rij, and R are, respectively, the covariant 
derivative, Ricci tensor, and Ricci scalar associated with 
the spatial metric 7^ . Finally, the trace of the extrinsic 
curvature is denoted K = K\. 

The task of constructing initial data for a Cauchy evo- 
lution via Einstein's equations requires that we decom- 
pose the constraints in such a way that we can specify 
how the constrained, gauge, and dynamical degrees of 
freedom are associated with the initial data. In this pa- 
per, we are primarily interested in initial data associated 
with systems in quasiequilibrium. Because of this, it is 
natural to use the the conformal thin-sandwich decom- 
position of the constraints |3, |3j. This decomposition is 



particularly useful in this situation because quasiequi- 
librium is a dynamical concept, and this decomposition 
retains a close connection to dynamics that is lost in most 
other decompositions of the constraints (cf Refs. 0,3). 

The conformal thin-sandwich decomposition employs a 
York-Lichnerowicz conformal decomposition of the met- 
ric and various other quantities 0, 0, Q- The conformal 
factor, ip, is defined via 



(7) 



where 7^ is a "conformal metric" . The time derivative 
of the conformal metric is introduced by the definitions 



7 tJ «ij = 0. 



(8) 
(9) 



From this, it follows that the tracefree extrinsic curvature 
A* ee K l i - WiK takes the form 



2a 



(10) 



where 



= ,1,-6 



u VJ = Ufcz7 ifc 7 j ' 



ip a is the conformal lapse function, and 



Furthermore, (LV) is the conformal- 
Killing (or longitudinal) operator acting on a vector, de- 
fined by 



2V ( ^-) 



(11) 



where is the covariant derivative compatible with 7^ . 

Notice that this decomposition of incorporates the 
kinematical variables of the 3+1 decomposition, that is, 
the shift vector j3 l and the lapse function a through the 
conformal lapse a. It also includes the trace-free time 
derivative of the conformal metric, Uy. Below, the con- 
formal tracefree extrinsic curvature will be useful, 



A ij ee ip 10 A lj = — 
v 2a 



(12) 



Within the conformal thin-sandwich formalism, one 
must specify 7^, Uij, K, and a. With these quantities 
defined, the Hamiltonian (JjjJ and momentum con- 
straints take the form of a coupled set of elliptic equations 
that determine ip and /?'. In terms of our conformally de- 
composed variables, the Hamiltonian constraint (JSJ) can 
be written 

V 2 V> - - j2^ 5 K 2 + |V -7 'Aij&i = 0, (13) 

where R is the Ricci scalar associated with 7^, and the 
momentum constraint © as 

Vj ~ - V, = 0. (14) 

The freely-specified data includes the conformal met- 
ric, jij, and its time derivative, «y = 9t7ij) a s well as 
the trace of the extrinsic curvature, K , and the confor- 
mal lapse, a. It is possible, and desirable, to make the 
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set of freely-specified data more symmetric by choosing 
to specify the time derivative of the extrinsic curvature 
instead of the conformal lapse. This is possible because 
these two quantities are related by the trace of Eq. 
The resulting equation is an elliptic equation for the con- 
formal lapse that is coupled to both the Hamiltonian and 
momentum constraints, l|13[) and Q14JI. This equation can 
be written as 

= -tA 5 (d t K - f3 k V k K^j . (15) 

The statement made earlier, that the conformal thin- 
sandwich decomposition has a close connection to dy- 
namics is now clear. Not only does this decomposition 
incorporate the kinematical variables of the 3+1 split- 
ting, but fully half of the freely-specifiable data consist of 
time derivatives of fundamental fields. In particular, we 
are free to choose the conformal metric and the trace of 
the extrinsic curvature (7^ and K) and the time deriva- 
tives of these fields (d t %j = Uij and d t K). 

Remaining to be determined are the conformal factor, 
tp, the conformal lapse, a, and the shift vector, fj 1 . These 
are determined by solving Eqs. l(T3|) . ((111) . and (|T3|l as 
a coupled set of elliptic equations. Formulating a well- 
posed system requires that we impose boundary condi- 
tion. Typically, these systems are solved under the as- 
sumption that the spacetime is asymptotically flat. If we 
let r denote a coordinate radius measured from the loca- 
tion of the center of energy of the system, then asr^oo 
we have that 

V>| = 1, (16a) 
^U^^oxr) 1 , (16b) 



£!o is the orbital angular velocity of a binary system, 
or the rotational angular velocity of a single object, as 
measured at infinity. The boundary condition on the 
shift is chosen so that the time coordinate, i M = an M +/3 M , 
is helical and tracks the rotation of the systemQ, H, 0- 
If we wish to consider systems with one or more black 
holes, and if we excise the interior of the black hole to 
avoid difficulties with singularities, then we will also need 
to impose boundary conditions on the excision surfaces. 
This is the topic of the next section. 

III. BLACK-HOLE EXCISION BOUNDARY 
CONDITIONS 

The physical content of initial data depends on the 
choices made for the initial-data decomposition scheme, 
freely specifiable data, and the boundary conditions. 
Therefore, it is important to choose boundary conditions 
that are motivated by, or at least compatible with, the 
sort of initial data that we wish to construct. 



The first attempts to impose boundary conditions on 
black hole excision boundaries were based on topologi- 
cal arguments 0, 0> IS E3- By demanding that the 
initial-data hypersurface consist of two identical (isomet- 
ric) asymptotically flat hypersurfaces connected together 
at a number of spherical excision surfaces (one for each 
black hole) , it is possible to show that the surfaces where 
the isometric sheets connect are fixed point sets of the 
isometry. This condition is enough to determine either 
Dirichlet or Neumann boundary conditions at the exci- 
sion surface for any fields that are present. 

Boundary conditions based on this idea have been used 
successfully for generating general black-hole initial data 
using various initial-data decompositions [Til ITEl Pill W \ . 
Their first use in conjunction with the conformal thin- 
sandwich decomposition |9l 1 1 8| was only partially success- 
ful due to an unavoidable constraint violation. The diffi- 
culties with this approach were outlined in Ref. 1] , where 
an alternative approach of using quasiequilibrium bound- 
ary conditions was first outlined. In this section, we will 
refine, and in subsequent sections test, this approach. 

In constructing initial data on a spacelike hypersur- 
face, we cannot have knowledge of the event horizon that 
is typically used to define the surface of a black hole. 
However, we can identify the apparent horizon of a black 
hole, defined as the outermost marginally outer-trapped 
surface. A marginally outer-trapped surface (MOTS), in 
turn, is a surface on which the expansion, 8, of the family 
of outgoing null rays, vanishes everywhere. 

In this paper, we are interested in the situation in 
which each black hole is in quasiequilibrium. The as- 
sumptions required to enforce this are essentially the 
same as those required of an "isolated horizon" (cf 0, 
HJIHl). To ensure that the black hole is in quasiequi- 
librium, we enforce the following conditions. First, we 
demand that the expansion 9, vanish on the excision sur- 
face, S, thus forcing the boundary to be an apparent 
horizon: 

0\ S = O. (17) 

Next, we require that the shear <7^„ of the outgoing null 
rays also vanish on the excision boundary, 

<V,| S = 0. (18) 

Consider the family of null geodesies threading the ap- 
parent horizon to the future of our initial-data slice that 
are tangent to k 11 on S. Raychaudhuri's equation for null 
congruences, 

C k 9 = -\6 2 - a^a^ + w^w"" - R^k v , (19) 

together with Eqs. (jl7(l and l|18(l are sufficient to imply 
that 

c k e\ s = o. (20) 

That is, initially, the apparent horizon will evolve along 
/c M . Note that ui^ is the twist of the congruence, which 
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vanishes because the congruence is surface forming. Also, 
we assume that there in no matter on S, so = 0. 

While conditions (fTTjl . flE) and $2® are coordinate 
independent, our next and final demand breaks precisely 
this coordinate freedom: We demand that the coordinate 
location of the apparent horizon does not move initially 
in an evolution of the initial data. 

As we show in the subsequent sections, the require- 
ments listed so far yield four conditions that can be im- 
posed on the initial data at the excision boundary. How- 
ever, there are five coupled initial-data equations that 
must be solved in the conformal thin-sandwich approach. 
When quasiequilibrium black-hole boundary conditions 
were first derived, a fifth condition was considered 0. In 
particular, the condition that Ck9 = was considered, 
where 9 is the expansion of a family of ingoing null rays, 
k^ . As we will show below, this fifth condition cannot be 
used as a boundary condition even though it is satisfied 
for a stationary black holel 

In the remainder of this section, we will derive bound- 
ary conditions for black hole excision surfaces based on 
the demands outlined above. A good portion of the fol- 
lowing derivation appeared previously 1]. However, be- 
cause of a change in notation, and more importantly in 
a few sign conventions, we include the full derivation be- 
low. 



where t/ M „ is the full spacetime metric. Because these ten- 
sors T,^ and T,^ are spatial, we will use spatial indices 
below. To simplify the definitions that follow, we will 
introduce various projections of Kij along and normal to 
the excision boundary S: 

Jij ee hthtKM, (26) 
Ji = h\s l K u , (27) 
J ee h ij Jij = h i] K tj . (28) 

We can then simplify Eqs. (|21fl and lt251) to 

Sij = -^(Jij-Hij) and = —j%(J%j +Hij). (29) 

Now, we define the expansion of outgoing null rays, 9, 
and ingoing null rays, 9, via 

9 EE ft«Sy =--^(J-H), (30) 

6 = = -±(J + H). (31) 

Finally, we define the shear of the outgoing null rays, <xy, 
and ingoing null rays, via 

o~ij = Ynj — 7^hij9, (32) 

<7jj = Y^ij (33) 



A. Geometry of the excision boundary 

We demand that the excision boundary surface, S, be 
a spacelike 2-surface with topology S 2 and define s l to be 
the outward pointing unit vector normal to the surface. 
In this case, we define outward with respect to the black 
hole (not the domain), so that s % points toward infinity. 
The 4-dimensional generalization of s l has components 
s M = [0, s 1 } obtained from the condition that s^n^ = 0. 

The metric, h%j } induced on S by 7^ is given by 

Kj EE Jij - SiSj. (21) 

We also define the extrinsic curvature, Hij, of S embed- 
ded in the 3-dimensional spatial hypersurface as 

Hy ee ^V (fe ^) = (22) 

Naturally associated with S are two sets of null vectors: 
a set of outgoing null rays, and a set of ingoing null 
rays, fc M , defined by 

k^= -^(n^ + s**) and U l = -±(n^ - s»). (23) 

Associated with each set of null rays is an extrinsic curva- 
ture of S as embedded in the full 4-dimensional manifold. 
These are defined as 

E pi/ ee ±h a XC k g af) , (24) 
£„„ ee \hlhiL k Q afs , (25) 



B. Quasiequilibrium boundary conditions 

With the definitions of Sec. lIII Al we can now evaluate 
the demands we made earlier in Sec. IIIII and translate 
them into boundary conditions for the conformal thin- 
sandwich equations. In order to express these as useful 
boundary conditions, we must write them in terms of 
the variables of the conformal thin-sandwich approach. 
We must also make connection with the global notion of 
quasiequilibrium that is closely associated with an ap- 
proximate helical Killing vector. 

A spacetime that is in true equilibrium is said to be 
stationary and has two Killing vectors of interest: a time- 
like Killing vector, d/dto, and a spatial Killing vector 
associated with rotational symmetry, d/d/fio- If fio de- 
notes the angular velocity of a spinning object or sys- 
tem as measured at infinity, then the linear combina- 
tion d/dto + £lod/dcf>o is referred to as the helical Killing 
vector. If a system, such as a binary, is in a state of 
quasiequilibrium, there arc in general no vector fields 
similar to d/dto or d/d(f>o that are even close to being 
Killing vectors. But there will be a helical vector field 
that is an approximate Killing vector of the spacetime. 
If we let this approximate Killing vector field define our 
time vector, i M , and thus our time coordinate t, then we 
will have d/dt ss for fields in this spacetime. Within 
the 3+1 decomposition, we write the time vector as 

= an fi +0". (34) 



5 



Our desire for t^ to represent an approximate helical 
Killing vector is, therefore, the reason for our condition 
on the shift at infinity, Eq. (|16bjl . 

Now we consider the demand that the apparent hori- 
zon should initially not move in an evolution of the quasi 
equilibrium initial data. Because of Eq. I(20|l . the appar- 
ent horizon initially coincides with the null surface gen- 
erated by k^. In order for the coordinates to track this 
null surface, the time- vector of the evolution, t^, must lie 
in this null surface. This requires that 



t hu 



0. 



(35) 



Substituting Eqs. and J23) into Eq. (J33J), and re- 
calling that the shift vector is spatial, fi^n^ — 0, yields 



(36) 



This equation is often referred to as the Killing-horizon 
condition. We split the shift vector into its component 
normal to the surface, (3 ± , and a vector tangent to the 
surface, /3/, defined by 



0± 



(37) 
(38) 



With these definitions, we see that Eq. I|36|) is a condition 
on the normal component of the shift, 



0± 



(39) 



The component of the shift tangential to the excision 
surface 5, ft/ remains unconstrained so far. This makes 
sense, because fixing the location of a surface does not 
restrict motion within this surface. We can gain insight 
into the relevance of /3y 1 by considering a stationary Kerr 
black hole. 

The Kerr spacetime has two Killing vectors of interest: 
A timelike Killing vector, d/dto, and a spatial Killing 
vector associated with rotational symmetry, d/dipo. The 
null generators of the horizon are given by 



_d_ 



n 



H 



d 



(40) 



where is the angular frequency of the horizon. If we 
introduce a fiducial helical Killing vector, 



d_ 



d 



for some fi, then, of course, on the horizon, 

d 



(41) 



(42) 



Now consider a hypersurface through Kerr to which 
<9/<90o is always tangent, e.g., the usual Kerr-Schild slice. 
To make the connection with the usual 3+1 decomposi- 
tion straightforward, we normalize k and £ by choosing 



their time components to be one (fc* = £* = 1). If we 
choose the vector £ as the time vector of an evolution, 
then the last term in Eq. 1)42(1 corresponds precisely to 
/3|| Z — this term is tangent both to the horizon and to 
the hypersurface. For the choice (3^ = 0, the last term in 
Eq. (|42|l would be absent, i.e. SI = £Ih- In this case the 
black hole is corotating with the coordinate system, and 
the generators of the horizon, fc M do not twist relative to 
the helical Killing vector, £ fJ -. Conversely, a non-rotating 
black hole (£Ih — 0) would have a tangential shift of 



(43) 



For a binary black hole in quasiequilibrium, neither 
d/dto nor d/dtfio exist as separate Killing vectors. How- 
ever, based on the discussion above, we expect that (3/ 
is still connected to the rotation of the black hole. In the 
case (3/ = 0, the horizon generators have no twist rel- 
ative to the helical Killing vector (which coincides with 
the time- vector) , corresponding to corotating black holes. 
Any rotation, (3 n l oc (d/d^s) 1 , where (d/dfisY lies in the 
surface S would impart additional rotation on the black 
hole. Below we will make these notions precise. 

Having obtained a boundary condition on (3 ± , we now 
turn our attention to Eqs. (|17fl and l|18f) . We need to 
consider the horizon boundary S in the conformal space. 
The conformal transformation on 7y J7J) induces a nat- 
ural conformal weighting for hij and for the unit normal 
to S, 



- 2 s\ 



(44) 
(45) 



If we also define D % as the covariant derivative compatible 
with hij , then without loss of generality, we can express 
the expansion of the outgoing null rays, 9, as 



V2 



h lJ V l s J ; +4§ fc V fe hn/>- if J) , 



(46) 



and the shear of the outgoing null rays, a lJ , as 



— <H lj - ±h i3 H) 11- — 



1 f 

V2 a 



-4 r 



(47) 



\{W k h? t u u - \~ti r h M u u ) 



It is now clear how to obtain the remaining bound- 
ary conditions. By applying condition lfT7l) to Eq. l(4^|) . 
we obtain a boundary condition that forces an excision 
boundary to be an apparent horizon (or MOTS). The 
condition is 



s^Vfeln^ 



i> 2 J 



(48) 
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and it takes the form of a nonlinear Robin-type boundary 
condition on the conformal factor, ip. Finally, if we recall 
the condition from Eq. (j39JI and that we have chosen 
u l i = everywhere, and apply condition 118fl to Eq. (|47|l . 
we obtain a condition that restricts the form of (3^ . The 
condition is that 



D"0 



= 0. 



(49) 



This shows that the components of the shift that are 
associated with the spin of the black hole must be pro- 
portional to a conformal Killing vector of the conformal 
metric, h%j, defined on the 2-dimensional excision surface. 

This condition is quite remarkable. Recall that any 
2-surface that is topologically S 2 is conformally equiva- 
lent to the unit 2-sphere. If you consider a unit 2-sphere 
embedded in a flat 3-dimensional Euclidean space, then 
there is a family of rotational Killing vectors, £*, associ- 
ated with any rotation axis passing through the center of 
the 2-sphere. Because ^ l fii = on the 2-sphere, where 
n l is the unit normal vector on the 2-sphere, we see that 

trivially form a family of 2-dimensional vectors tan- 
gent to the 2-sphere and that these are Killing vectors 
of the metric on the unit 2-sphere. But, the Killing vec- 
tors associated with any metric are also conformal Killing 
vectors of any metric conformally related to it. So, if tp 
represents a conformal transformation such that if hij is 
the metric of the unit 2-sphere, then will satisfy the 



conformal Killing equation on hij. Thus, 

/V = fW 4 



(50) 



will satisfy Eq. I|49|) . with fl r being an arbitrary param- 
eter. The freedom left in Eq. i|5U|) is precisely what is 
required to parameterize an arbitrary spin on the black 
hole. The parameter Q r is associated with the magni- 
tude of the rotation or spin of the black hole, whereas 
the axis of rotation of £ z is related to the orientation of 
the spin. Of course, Q r does not correspond directly to 
the rotational angular velocity of the black hole. From 
the discussion leading to Eq. (|43|l . it is clear that O r = 
corresponds to a black hole that is corotating with an ap- 
proximate helical Killing vector t^ and thus represents a 
black hole that is rotating as seen by an inertial observer 
at infinity. In order to construct a black hole that is not 
rotating as seen from infinity, we need to choose a shift 
that is similar in form to Eq. (|43|l . It seems reasonable to 
choose Sl r = fio and to pick the conformal Killing vector, 
£ l from the family that corresponds to rotation about an 
axis that is perpendicular to the plane of the orbit. 

To summarize, the quasiequilibrium conditions defined 
in Eqs. H17|) and l|18l) define boundary conditions on the 
conformal factor, ip, via Eq. I|48l) and on the shift vector, 
P\ via Eqs. $7$, (gHJl, ®, and (JSDJ. These total to four 
of the five necessary boundary conditions for solving the 
coupled elliptic equations associated with the conformal 
thin-sandwich equations. Missing is a condition on the 
conformal lapse, a. 



C. Boundary conditions on the lapse function 

A possible boundary condition on the lapse was de- 
rived in Ref. pj. The condition was essentially based on 
the reasonable quasiequilibrium condition that CuO = 0. 
Notice that this is the change in the expansion of ingoing 
null rays, fc M , with respect to the outgoing null congru- 
ence. The resulting boundary condition takes the form 



&BC 



0. 



(51) 



where £ is a nonlinear operator that is elliptic within the 
surface S (see Eq. (84) of Ref. Q for a precise descrip- 
tion) . This condition is satisfied on the horizon of a sta- 
tionary black hole and it seemed to supply a reasonable 
boundary condition for the lapse to be used in conjunc- 
tion with the previously defined boundary conditions for 
the conformal factor and shift vector. 

We implemented the full set of boundary conditions 
within the code described in Ref. (2^ and attempted to 
solve the full set of conformal thin-sandwich equations 
for the case of a single nonrotating black hole. It became 
clear immediately that the iterative solutions would not 
converge in general. Interestingly, if the analytic solution 
for an isolated black hole was supplied for the starting 
point of the iterations, then the equations and boundary 
conditions were satisfied to truncation error, but the iter- 
ations were at best only marginally stable. Furthermore, 
if any one of the boundary conditions on tp, f3 ± , or a were 
replaced by an arbitrary Dirichlet or Neumann boundary 
condition, then the iterative solution was convergent to a 
solution representing a static black hole and the omitted 
boundary condition was satisfied. 

This clearly indicates that the set of boundary condi- 
tions including Eq. (|5T|) is degenerate and leads to an 
ill-posed elliptic system. We can understand the na- 
ture of the degeneracy by considering the family of time- 
independent maximal slicings of Schwarzschild[23l I24L 
|25| . The line element for the spatial metric, lapse, and 
shift vector are 



ds 2 



= 1- 



2M C 2 



R 



R 4 



2M CP 
~R~ + R^' 



0« = 



C 2M CP 

R 2 \ li + R i ' 



dR 2 + R 2 dfl 2 , (52a) 
(52b) 
(52c) 



where R is the usual Schwarzschild "areal" radial coordi- 
nate, M is the mass of the black hole, and C is a constant 
parametrizing the family of maximal slicings. In spher- 
ical coordinates, the extrinsic curvature takes the form 



C_ 



-2 
1 
1 



(52d) 
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For < C/M 2 < « 1.299 the maximal spatial slice 

extends from spatial infinity, through the black-hole inte- 
rior, and to the second spatial infinity of the maximally 
extended Schwarzschild geometry. When C/M 2 = 0, 
we recover the standard Schwarzschild maximal slice 
that passes through the bifurcation point. And, for 

C/M 2 > yf|, the maximal spatial slice extends from 

spatial infinity, through the black-hole horizon, and ends 
on the singularity. 

If desired, the family of maximal slicings of 
Schwarzschild can be be rewritten in terms of an isotropic 
radial coordinate, r. It is then easy to verify that the 
boundary conditions (gSJ, (|2HJ>, and (|5T1) are satisfied 
on the horizon for any value of C/M 2 . So, we see that 
while the set of boundary conditions proposed in Ref. p) 
hold for a time-independent configuration, they do not 
uniquely fix the spatial slicing. It is in this way that 
they are degenerate. 

If we consider the value of the lapse on the horizon, we 
find that 



a(r n ) 



C 
M 2 



(53) 



Where rn denotes the location of the horizon in isotropic 
coordinates. Note that r n ^ M /2 unless C/M 2 = 0. If, 
instead of using the lapse boundary condition l]5ip. we 
simply fix a Dirichlet value for the lapse, then we find 
that we have effectively chosen a value of C /M 2 and thus 
a particular maximal slicing of Schwarzschild. Similarly, 
we find that the family of mixed boundary conditions 



da 
Or 



= A« 

r 



(54) 



corresponds to a slicing choice of C/M 2 = \/ A 2 + 4 — A, 
where A is any real number. It is clear that any reason- 
able choice of a Dirichlet, Neumann, or mixed boundary 
condition on the lapse will uniquely fix a particular max- 
imal slicing of Schwarzschild and effectively break the 
degeneracy. 

As mentioned previously, our numerical investigations 
have shown that we could also have chosen to use the 
lapse boundary condition (|5T)l and instead fix either tjj 
or (3 ± on the horizon via a Dirichlet, Neumann, or mixed 
boundary condition. One reason to choose to replace the 
lapse condition l|51|l. as opposed to the boundary con- 
ditions on the conformal factor l|48|) or the shift l|39l) . is 
that the lapse boundary condition is much more complex. 
However, there is a more fundamental reason to choose 
to replace the lapse condition. The degeneracy that we 
must eliminate is in the choice of the spatial slice which 
is a choice of the initial temporal gauge. The lapse func- 
tion fixes the evolution of the temporal gauge. Therefore 
it seems reasonable that we should consider the choice of 
the lapse boundary condition as part of the initial tem- 
poral gauge choice. It is customary to view the choice of 
the trace of the extrinsic curvature as fixing the initial 



temporal gauge. However, this is apparently not suffi- 
cient within the conformal thin-sandwich approach when 
interior boundaries are present. 

This last assertion, that the lapse boundary condition 
must be chosen as part of the initial temporal gauge 
choice, is supported by the behavior of a very special so- 
lution of Einstein's equations - the static Schwarzschild 
solution. However the same generic behavior is seen in a 
broad class of examples as we will outline below. 



IV. QUASIEQUILIBRIUM SOLUTIONS FOR A 
SINGLE BLACK HOLE 



As we have seen in Sec.[Hl the conformal thin sandwich 
equations require specification of free data, which are the 
conformal metric jij and its time derivative Uij , as well as 
the mean curvature K and its time derivative d t K. More- 
over, boundary conditions are required on the variables 
being solved for, the conformal factor ip, the shift f3 l and 
the lapse a. The quasiequilibrium approximation fixes a 
good portion of these choices, namely Uij = dtK = 0, as 
well as the the following boundary conditions at the ex- 
cised regions: The apparent horizon condition Eq. (|48|l on 
tp, the null horizon condition Eq. H39fl on ±) and Eq. H5U|) 
on /3||\ At spatial infinity, the boundary conditions are 
straightforward and are given by Eqs. (|16(l . 

Before the conformal thin-sandwich equations can be 
solved, we have to choose the remaining quantities which 
are not fixed by the quasiequilibrium framework. These 
are jij, K, and an inner boundary condition on the lapse 
a. Furthermore, we have to choose the shape of the ex- 
cised regions, S. As argued above, and as we confirm 
below, the lapse boundary condition and K are part of 
the temporal gauge choice. It thus remains to choose 7y 
and the shape of the excised regions, S. 

In this work, we will assume that the conformal three- 
geometry is fiat, and we will always excise exact spheres. 
These choices are not motivated by physical considera- 
tions, and they will affect the quality of the quasiequi- 
librium solutions we obtain in this paper. For exam- 
ple, the Kerr spacetime does not admit conformally flat 
slices |2(||27J. Therefore, when we solve for a rotating 
black hole, our initial-data sets will not exactly represent 
a Kerr black hole, but will rather correspond to a per- 
turbed Kerr black hole, which will settle down to Kerr. 
We stress that this failure of the initial-data sets con- 
structed here to represent Kerr is not inherent in the 
quasiequilibrium method, but is caused by our choices 
for 7^ and S. With the appropriate choices for 7^ and 
iS, the quasiequilibrium method can reproduce exactly 
any time-independent solution of Einstcins equations. In- 
deed, for single black holes, better choices for 7^ and S 
are easily obtained from stationary analytic solutions of 
Einstein's equations, for example based on Kerr-Schild 
coordinates. While such a choice certainly leads to single 
black hole initial-data sets closer to the true Kerr met- 
ric, it is not clear how to generalize to binary black hole 
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configurations. A widely used approach superposes sin- 
gle black hole quantities to construct binary black hole 
initial data (e.g., [3 S El ED, EH)- However, be- 
cause of the nonlinear nature of Einsteins equations, and 
since often, the black holes are separated by only a few 
Schwarzschild radii, the superposition introduces uncer- 
tainties that may be large |23 and that have not yet been 
adequately quantified. In this paper, rather than using 
superposition, we start with choices for 7y and S that are 
not optimal for single black hole spacetimes, but that are 
equally well suited for binary black hole configurations. 
We then use the single black hole solutions to quantify 
the effects of our approximation. 

We solve the conformal thin-sandwich equations with 
the pseudo-spectral collocation method described in [3^] . 
For the single black hole spacetimes, typically, two spher- 
ical shells arc employed. The inner one ranges from the 
excised sphere to a radius of ~ 20 and distributes grid- 
points exponentially in radius. The outer shell has an 
outer radius of typically ~ 10 10 , and employs an inverse 
mapping in radius, which is well adapted to the 1/r— 
falloff of many quantities. The fifth elliptic equation (for 
the lapse function) is coded in the form of Eq. 115|) . i.e. 
as an equation for ip 7 a — aip. Therefore, we formulate 
the lapse boundary condition as a condition on aeip. Fi- 
nally, we note that we always solve the three-dimensional 
initial value equations, even in cases which have spher- 
ical or cylindrical symmetries like the single black hole 
solutions. 



A. Spherical symmetry 

We begin by solving for spherically symmetric initial- 
data sets that contain one black hole. In spherical sym- 
metry, the assumption of conformal flatness is no restric- 
tion, because any spherically symmetric metric can be 
made conformally flat through an appropriate radial co- 
ordinate transformation. For example, a hypersurface 
through the Schwarzschild spacetime of constant Kerr- 
Schild coordinate time has the induced metric 



ds 2 = ( 1 + ) dR 2 + R 2 dtt 2 , 



(55) 



where R denotes the areal radius. Here, the coordinate 
transformation p| 



R 



brings the induced metric into conformally flat form, 

ds 2 = V4 (r 2 dr 2 + rW) , (57) 



with ijj KS = \/R/r. 

In order to support our claim that K and the lapse 
boundary condition merely represent a coordinate choice, 



we solve the quasiequilibrium equations for several, es- 
sentially arbitrary, choices for these quantities. For the 
mean curvature, we choose 



K = 0, 
2M 



K 



1 



2M 
~R 



-3/2 



1 



(58a) 
(58b) 

~R~) ' ^ 



where in the last case, R is given implicitly by Eq. 15611 . 
Equation l|58c[) represents the mean curvature for a Kerr- 
Schild slice of the Schwarzschild spacetime with mass M. 

For the lapse boundary condition at the excised 
spheres, we use 



daip 



dr 
daip 



dr 



am 



0. 

cap 
"27 
1 

2' 



(59a) 

(59b) 

(59c) 
(59d) 



The last condition, Eq. I|59d|l is correct for the Kerr- 
Schild slice. 

We now compute twelve initial-data sets, combining 
any of the choices for K with any of the lapse boundary 
conditions. In order to fully recover the Kerr-Schild slice 
for the choices (|58cll and (|59djl , we set the radius of the 
excised sphere to 



I R=2M 



1 + V2 



„2-2%/2 



M. 



(60) 



For each of the initial-data sets, we compute the 
residual in the Hamiltonian and momentum constraints, 
Eqs. © and ©. We compute ADM quantities of the 
initial-data set by the standard integrals at infinity in 
Cartesian coordinates, 



^ADM — 



J (?) = ~ 



1 

16tt 
1 

8^r 



(61) 
(62) 



For the ^-component of the linear ADM-momentum, 
£ = e x in Eq. (|62|l . The choice £ = xe y — ye x yields 
the z-component of the ADM-like angular momentum as 
defined by York [3j| • We also compute the irreducible 

mass 



A A 



H 



16tt 



(63) 



where we have approximated the area of the (unknown) 
event horizon by the area Aah of the apparent horizon, 
and the Komar mass, 



MK 47T ./_„ 



(%a-l3 k K lk ) d 2 S 3 



(64) 
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FIG. 1: Solution of the quasiequilibrium equations recovering 
the Kerr-Schild slicing of Schwarzschild. Plotted are the max- 
imum values of Hamiltonian and momentum constraints and 
of time-derivatives, as well as the deviation of Mi rr , E A um and 
Mk from the analytical answer M. N r is the radial number 
of collocation points in each of the two spherical shells. 



Finally, we evaluate the time-derivative oiKij by Eq. |@}, 
evaluate the residual cxbc of the quasiequilibrium lapse 
condition, Eq. (|51|l . and compute 

d t \nip = ~(aK (65) 

which follows from the trace of Eq. . 

Figure ^ presents a convergence plot for one of the 
twelve cases, Eqs. (I58cf) and (I59d|l . This case recovers 
the usual Kerr-Schild slice with mass M . The resid- 
ual of the Hamiltonian and momentum constraints de- 
crease exponentially with resolution, and the three dif- 
ferent masses M; rr , E Aum and M K all converge to the ex- 
pected result, M . Furthermore, the time-derivatives ex- 
ponentially converge to zero, and ADM linear and angu- 
lar momenta converge to zero, too. The vanishing time- 
derivatives indicate that the quasiequilibrium method 
constructs lapse and shift along the timelike Killing vec- 
tor of the Schwarzschild spacetime. 

For different choices of K or for the lapse BC, we find 
that the masses are no longer exactly unity. However, for 
all choices of K and lapse BC, we find to within trunca- 
tion error, that the three masses agree, 

M ilT = E ADM = M K , (66) 

and that all time-derivatives (and the lapse condition 
Eq. jSJ) vanish: 

dtjt = dtKij = a BC = 0. (67) 



TABLE I: Spherically symmetric quasiequilibrium initial-data 
sets. Given are the irreducible mass, ADM-energy and Komar 
mass (these three quantities are found to be identical to within 
truncation error). The last column gives an upper bound 
on the deviation from zero of Hamiltonian and momentum 
constraints, time-derivatives of tp and Kij, as well as the lapse 
condition Qsc 
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These findings are summarized in Table [I] These runs 
indicate that any (reasonable) choice for the mean cur- 
vature K and the lapse BC recovers a slice through 
Schwarzschild with time- vector along the timelikc Killing 
vector. 

For the maximal slices, K = 0, the different lapse 
BCs choose different parameters C/M 2 in the family of 
maximal slicings, Eqs. I|52a|l - lj52d|l . The boundary con- 
ditions in Eqs. I|59a() - (|59d|l correspond, respectively, to 
C/M 2 = 2/3(^13 - 1), C/M 2 = 4/3, C/M 2 sa 1.2393, 
and C/M 2 w 2.4905. Based on the results of Table □ 
we conjecture that for any (reasonable) function K(r), 
there exists a one-parameter family of spherically sym- 
metric slicings, which extend from the horizon to spatial 
infinity. 

In situations with less symmetry like binary black 
holes, we prefer Neumann or Robin boundary conditions 
on the lapse (Eqs. |59aj and |59b| L because they allow 
the lapse on the horizon to respond to tidal deformations. 
Furthermore, as Table |U confirms that the choice of mean 
curvature plays a marginal role, we will concentrate on 
the most obvious choice, maximal slicing K = 0, below. 

In summary, the quasiequilibrium method is singularly 
successful for spherically symmetric spacetimes: There 
exists a natural choice for 7^ (the flat metric), that, to- 
gether with any (reasonable) choices for K and the lapse- 
boundary condition yields a slice though Schwarzschild 
and the timelike Killing vector that results in a com- 
pletely time-independent evolution. 

B. Single rotating black holes 

In Sec. IIIII we explained that the tangential component 
of the shift-vector on the excised surface induces a rota- 
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FIG. 2: Spinning single black hole initial-data sets. 

tion on the hole, cf. Eq. I|50|) . We now test this assertion 
by constructing initial-data sets for single rotating black 
holes. 

We set the mean curvature and the lapse BC by 
Eqs. (|58a|l and l|59a[l . and we continue to use the con- 
formal flatness approximation, 7^ = /y, and choose the 
excised region to be a coordinate sphere centered on the 
origin with radius r exc = 0.8594997. This value ensures 
a unit- mass black hole in the limit of no rotation for the 
choice of K and lapse BC. The shift boundary conditions 
encode the rotation. At infinity, we set j3 l — 0; on the 
horizon, Eq. (|50|l implies 



nix k e ijk , 



(68) 



where x k is the Cartesian coordinate separation of points 
on S to the center of the excised sphere. We choose 
parallel to the z-axis, and solve the initial value equations 
for different magnitudes of il r . For each solution, we 
compute the diagnostics mentioned in Sec. II V Al Figure[21 
presents the ADM-energy, irreducible mass and angular 
momentum of the obtained data sets. We see that, for 
small Q r , the angular momentum increases linearly with 
f2 r , as expected. 

As discussed early in Sec. lIVI our assumption of confor- 
mal flatness will necessarily introduce some errors when 
solving for a rotating black hole, because the Kerr met- 
ric does not admit conformally flat slices|26l l27j. Very 
interesting are therefore measures of the deviation of the 
quasiequilibrium initial-data sets to a slice through the 
exact Kerr spacetime. 

One such quantity is the maximum amount of energy 
that can potentially be radiated to infinity, 



E ra d — V^adm -PIe 



7 2 

_ ADM 

4ML 



(69) 



FIG. 3: Single rotating black hole initial-data sets: Deviations 
from the exact Kerr metric. 



For a stationary spacetime, E rac [ = 0. Another interest- 
ing question is how closely O r of Eq. (fFJHf corresponds to 
the angular frequency of the horizon. For a Kerr black 
hole with angular momentum J ADM and total mass £ ADM , 
the angular frequency of the horizon is given by |34| 



^ADM/^ADk 



2 + 2^/l-(J ADM / J E ADM ) 2 ' 



so that 



Ail = E ABM (Q r — Qh) 



(70) 



(71) 



measures the deviation of tt r from the angular frequency 
of the horizon. Figure presents these quantities. The 
maximum radiation content, E rac i, is proportional to 
(Mj rr f2 r ) 4 . For the binary black hole data sets we con- 
struct below in Sec. [3 the relevant angular frequency 
is the orbital angular frequency. We find below, that at 
the innermost stable circular orbit, M^Qq ~ 0.11. From 
Fig. |2| we find for this angular frequency, E ra d/ E AriM ss 
2 • 10~ 4 . This indicates that, when conformal flatness 
is assumed, we should not expect the fractional error 
induced in E ADM for systems with non-vanishing angu- 
lar momentum to be larger than ~ 10~ 3 . Furthermore, 
at Mi„il r = 0.11, Af2 pa 0.0025. From this we expect 
that the rotational state of each black hole, within the 
binary black hole configurations below, should deviate 
by at most one or two percent from the intended val- 
ues for corotating or irrotational holes. In Fig. [3] we also 
plot some time-derivatives assuming the initial data are 
evolved with the constructed gauge a,(3 l . These time- 
derivatives are proportional to (Mj rr f2 r ) 2 . However, their 
interpretation is more difficult, due to their gauge depen- 
dence and the difficulty of finding a meaningful normal- 
ization. 
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FIG. 4: Single boosted black hole initial-data sets 



C. Single boosted black holes 

A boosted single black hole in a comoving coordinate 
system appears time-independent. A well-known exam- 
ple is the boosted Kerr-Schild form of a Kerr black hole. 
In such comoving coordinates, the shift does not vanish 
at infinity, but approaches the boost velocity of the black 
hole, 

f3 l \ = v*. (72) 

We apply now the quasiequilibrium formalism to con- 
struct boosted black holes in comoving coordinates by 
using Eq. (|72[) as the boundary condition on the shift 
at the outer boundary. At the excised sphere, we set 
/3|| l | 5 = 0. Furthermore, we assume again conformal 
flatness, use Eqs. (|58a|l and l|59a[) to fix the mean cur- 
vature and the lapse BC, and excise a coordinate sphere 
with radius r exc . The remaining free parameter is the 
magnitude of the boost velocity, v. Figure 21 presents 
the ADM-energy, irreducible mass and P A dm/E A dm as a 
function of the boost-velocity. P A um/ E ADM is linear in v 
for small v, as it should be. As v approaches unity, i? ADM 
strongly increases. Figure |5] presents measures of how 
faithfully these initial-data sets represent a boosted sta- 
tionary black hole. The maximum radiation content E ra d 
grows as u 4 . At the ISCO we can estimate that v ~ 0.4 
where we find E ra d ~ 10~ 3 E ADM . In order to measure 
how well the special relativistic relation v = P/E is sat- 
isfied, we define 

P 

A 1 ADM /TO\ 

Av = v- ; (73) 

^ADM 

we find that Av oc v 3 for small v, with At; w 0.01 for 
v = 0.4. 
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FIG. 5: Single boosted black hole initial-data sets: Deviations 
from exact time-independence. 



V. QUASI-CIRCULAR ORBITS FOR 
BLACK-HOLE BINARIES 

In the case of a single black hole, if appropriate choices 
for the freely specifiable data and boundary conditions 
are made, then an exact equilibrium solution of the 
initial-data equations can be found. However for bi- 
nary black-hole configurations, no such true equilibrium 
or stationary state exists. This is a much more strin- 
gent test of the quasiequilibrium boundary conditions. 
In this section, we will examine the solutions for the case 
of equal-mass black-hole binaries that are either corotat- 
ing or irrotational. 

We will consider binary configurations over a range of 
separations. A black hole is represented in the coordinate 
system by an excised 2-surface. The relative coordinate 
sizes of various excised surfaces parameterize the relative 
sizes of the resulting physical black holes. The coordi- 
nate separations of the holes parameterize their physi- 
cal separation. These coordinate sizes and separations 
are measured in the coordinates associated with the cho- 
sen conformal metric. In this work, all excised surfaces 
are the surfaces of coordinate spheres. Furthermore, for 
the simple cases of corotating and irrotational binaries, 
equal-mass black holes are obtained by choosing excision 
surfaces for the two holes that have equal radii. 

Before solving the initial-data equations, we must 
make choices for the freely specifiable data. In all cases, 
we will make use of the quasiequilibrium assumptions 
on the free data that Uij — and dtK — 0. We also 
continue to use the approximation that the conformal 
three-geometry is flat (ie. that 7^ is a flat metric). The 
remaining free data is the trace of the extrinsic curva- 
ture, K . For this, we will consider two choices: maximal 
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slicing with K = 0, and a non-maximal slicing based on 
Eddington-Finkelstein slicing . 

In addition to the freely specifiable data, we must also 
fix the boundary conditions on the excision surfaces that 
correspond to the surface of each black hole and at the 
outer boundary of the computational domain. The outer 
boundary conditions were discussed at the end of Sec.lTTI 
The boundary condition on the shift, Eq. I|16b|) . contains 
a free parameter, Oo, that determines the orbital angu- 
lar velocity of the system. The value of this parameter is 
chosen by demanding that the ADM and Komar masses 
of the system must be equal 0,0,0] This is a quasiequi- 
librium condition that is satisfied by a single value of flo 
and places the binary in a nearly circular orbit. 

For the excision boundaries, we will use the apparent 
horizon condition given by Eq. I|48|l as a boundary condi- 
tion of the conformal factor ijj, and we will use Eq. I(39|l to 
fix the component of the shift that is normal to the exci- 
sion surface. Boundary conditions on the components of 
the shift that are tangent to the excision surface depend 
on our choice for the spins of the black holes. For the case 
of corotation, we will demand that (3/ — 0. However, the 
irrotational case requires a somewhat more complicated 
choice. 

The condition of quasiequilibrium requires that we 
choose the tangential components of the shift so that they 
have the form given in Eq. 1)50(1. For irrotational black 
holes in a binary, it is reasonable to choose the conformal 
Killing vector £ l so that it represents rotation about an 
axis that is orthogonal to the plane of the orbit. If we let 
n i r represent an angular velocity vector that is orthog- 
onal to the plane of the orbit, and if we use Cartesian 
coordinates for our flat conformal metric, then Eq. 15U|) 
can be written as 

/V = n?4e« fc , (74) 

where x± = x % — C % ± and C± is the Cartesian coordinate 
location of the center of either of two excision spheres. Fi- 
nally, we take the magnitude of £L l r to be equal to the or- 
bital angular velocity of the binary system as measured at 
infinity, Hq. We note that there is no rigorous proof that 
these choices lead to an irrotational binary system. How- 
ever, as argued in Sec. IIII Bl especially the paragraphs 
leading to Eq. (|43|) . these choices seem reasonable. 

Finally, we must choose a boundary condition on the 
lapse at the excision boundaries. For all choices of the 
black-hole spins and choices for K, we repeat the com- 
putations for three different lapse boundary conditions, 
namely Eqs. (|59a|l — (|59c|) . 

The conformal thin-sandwich equations are solved with 
the pseudo-spectral collocation method described in • 
The computational domain consists of one inner spherical 
shell around each excised sphere, which overlap 43 rect- 
angular blocks, which in turn overlap an outer spherical 
shell extending to r out = 10 9 . Fig. [5] shows convergence 
of this solver with spatial resolution for one typical con- 
figuration (separation d = 9, K = 0, corotating black 
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FIG. 6: Convergence of the elliptic solver for the binary black 
hole configurations. Plotted are the constraint violations in 
Hamiltonian and momentum constraint, and differences to 
the highest resolution solve. N is the cube-root of the total 
number of grid-points. 

holes). The calculations below are performed at a reso- 
lution comparable to N — 60, so that the discretization 
errors in E ADM and M K should be about 10 -6 . 



A. Maximal slicing 

1. Corotating binary systems 

We now compute initial-data sets corresponding to a 
binary black hole system in a quasi-circular orbit for 
many different separations. Figure [7| shows the bind- 
ing energy Ej, of the binary system as a function of the 
total angular momentum of the system for all the lapse 
boundary conditions. The binding energy is defined as 
Eb = £ AD m — m where E ADM is the total (ADM) energy 
of the system and m = mi + mi is the total mass of 
the system. For the quasiequilibrium numerical results 

described in this paper, we take tnwp, = \J as the 
irreducible mass of each individual black hole and A^ 2 
are the areas of the apparent horizon of each hole. The 
reduced mass is defined as a = mi " 12 . 

' m 

We note that the choice of the lapse boundary condi- 
tion has very little effect on the solutions. This is consis- 
tent with out assertion that the choice of the lapse bound- 
ary condition is part of the initial temporal gauge choice. 
The inset in Fig. \7\ shows a magnified view of the region 
where the black holes are closest to each other. Even 
in this region, the result of the different lapse boundary 
conditions are nearly indistinguishable. Because of this, 
subsequent plots displaying corotating maximal slicing 
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FIG. 7: Constant Mi rr sequence of corotating equal-mass 
black holes. Maximal slicing is used in these cases, and three 
different excision boundary conditions for the lapse were used. 



FIG. 8: Constant M- m sequence of corotating equal-mass 
black holes. Comparison of post-Newtonian EOB sequences 
with numerical maximal slicing results from HKV and this 
paper. 



results will only display one of these sequences. 

Figure [S] shows a comparison of the same data to anal- 
ogous results obtained by Grandclement et al.jljl] (la- 
beled CO:HKV-GGB in figure legends), and effective 
one-body post-Newtonian results as reported in Ref. [35[ . 
First, second, and third post-Newtonian (PN) results 
are displayed (labeled CO:EOB-lPN, CO:EOB-2PN, and 
CO:EOB-3PN respectively in figure legends). The 3PN 
results correspond to the approach labeled "3PN corot. 
A(u, a 2 )" in Table I of Ref. @. 

There is good agreement between all of the results at 
large separation. Also, it appears that the PN results 
are converging toward the quasiequilibrium numerical re- 
sults, even when the black holes are quite close to each 
other as seen in the figure's inset. We also see that the 
numerical results obtained by Grandclement et al.|l8| 
(hereafter GGB) differ only slightly from the quasiequi- 
librium results. As discussed in Ref. Q], the numerical 
solutions obtained by GGB must violate the constraints. 
The agreement seen in Fig. [S] lends support to the belief 
that the violation of the constraints is, in some sense, 
small and has a small impact on the physical content of 
the data. 

Another method for comparing data for the circular 
orbits of compact binaries is to examine the location of 
the inner-most stable circular orbit (ISCO). The 1SCO is 
not a well defined concept in general. However, in situ- 
ations where the dissipative effects of radiation reaction 
have been eliminated, an ISCO becomes more meaning- 
ful. The ISCO is defined in terms of a minimum of some 
appropriate energy. For corotating binary systems, true 
stationary configurations can exist (although they con- 
tain an infinite amount of energy in the form of gravi- 



tational radiation, see Ref. 36]). In this case, the min- 
imum of the total energy can be rigorous ly a ssociated 
with the onset of a secular instability [37], |38|. In the 
absence of true stationary configurations, this "turning 
point" method is still used to define the ISCO. 

In order to locate a turning point, one must have a se- 
quence of binary configurations spanning a range of sepa- 
rations. How this sequence is constructed is not uniquely 
defined. The ambiguity arises because of the lack of a 
fixed fundamental length scale in the problem. In the first 
work to construct sequences of black-hole binary initial- 
data sets representing circular orbit s[39|. the total mass 
of the black hole (defined in terms of the Christodoulou 
mass formula |40|l was used to normalize the sequences. 

GGB suggest another approach based on the thermo- 
dynamic identity |3^| 

d-E A DM = f^o dJ A DM- (75) 

This identity should be satisfied by a true stationary se- 
quence of corotating black holes. Let s denote some pa- 
rameter along a sequence of initial-data sets, and let e(s), 
j(s), and lu(s) denote the numerical values for dimen- 
sionless versions of the total energy, total angular mo- 
mentum, and orbital angular velocity at location s along 
the sequence. We are free to define a fundamental length 
scale x(s) along the sequence in any way we like, so long 
as we define the dimensionful total energy £ ADU (s), total 
angular momentum J A dm(s), and orbital angular velocity 
fio(s) consistently via, 

E ADM (s) = x(s)e(s), (76) 
JaomOO = X%s)3(s), (77) 
no(s) = x^sMs). (78) 
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Enforcing the identity (|75(l is sufficient to determine the 
change in x(s) between two points on the sequence. If we 
integrate along the sequence from a point s\ to another 
point S2, then we find that 

f f e'(s) -cj(s)j'(s) , I 
I i sa e(s)-w(a)j(s) J 

where a prime denotes differentiation along the sequence. 




FIG. 9: Plot of the irreducible mass of one black hole when the 
sequence of solutions is normalized to maintain AE = QdJ. 
Three different choices for the lapse boundary condition are 
shown for both corotating (upper plot) and irrotational (lower 
plot) black holes in an equal-mass binary. The vertical dashed 
line shows the approximate location of the ISCO defined by 
a minimum in Eh- The vertical dotted line shows the approx- 
imate location of the ISCO defined by a minimum in _Z?adm- 
In the irrotational case, this latter ISCO line is off the plot to 
the right. 

If the sequence is normalized via Eq. J75J, then the 
irreducible mass of one black hole, M- m = \m, is not 
necessarily constant along the sequence. The top half 
of Fig. shows Mj rr for a corotating quasiequilibrium 
equal-mass binary as a function of the orbital angular 
velocity. The length scale has been normalized so that 
M m = 1/2 at infinite separation. We confirm the finding 
of GGB that M; rr is nearly constant along the sequence. 
While there is a clear increase in the mass as the sep- 
aration decreases, this increase is small and appears to 
be of roughly the same order of magnitude as the differ- 
ences due to using different lapse boundary conditions. 
As we will see later, this behavior is not mirrored in the 
irrotational data. 

Yet another approach for normalizing the sequence is 
to demand that the individual irreducible masses associ- 
ated with the apparent horizons remain constant. This 
is a particularly convenient normalization for numerical 



work since it relies on a well defined and easily measured 
geometric quantity. Each of these choices for normaliz- 
ing an initial data sequence can effect the location of the 
ISCO, so it is important to use a consistent definition 
when comparing data. 

It is also important to clearly define which energy is 
being extremised when using a turning-point method to 
locate the ISCO. We can consider using the minimum in 
either the ADM energy E ADM or the binding energy Eh 
along any sequence to define the ISCO. From the defini- 
tion of the binding energy as Eh = E ADM — to, we see that 
the minima will not necessarily agree if m varies along 
the sequence. 

For the PN sequences, the ISCO is defined as the mini- 
mum in the binding energy E\, along sequences where the 
irreducible masses of the black holes remain fixed. For 
the PN sequences, Eq. (|75fl is identically satisfied as well, 
so this is equivalent to finding the minimum in -E A dm- 
This is not true for the quasiequilibrium numerical data. 

TABLE II: Parameters of the ISCO configuration for corotat- 
ing equal-mass black holes computed with the maximal slicing 
condition. Results are given for three different choices of the 
lapse boundary condition and two choices for the definition of 
the location of the ISCO. For comparison, the lower part of 
the table lists results of Refs. [3 HI El; "PN standard" 
represents a post-Newtonian expansion in the standard form 
without use of the EOB-technique. 



Lapse BC ISCO min. 


mflo Eb/m J/m 2 


d(^) = ADM 

a 7- 

E b 


0.105 -0.0165 0.844 
0.107 -0.0165 0.844 


aip = ± ADM 
E b 


0.106 -0.0165 0.843 
0.107 -0.0165 0.843 


d(^, = ADM 

ar Zr 

E b 


0.106 -0.0165 0.843 
0.107 -0.0165 0.843 


HKV-GGB 


0.103 -0.017 0.839 


IPN EOB 
2PN EOB 
3PN EOB 


0.0667 -0.0133 0.907 
0.0715 -0.0138 0.893 
0.0979 -0.0157 0.860 


IPN standard 
2PN standard 
3PN standard 


0.5224 -0.0405 0.621 
0.0809 -0.0145 0.882 
0.0915 -0.0153 0.867 



Table IITI displays the dimensionless orbital angular ve- 
locity, binding energy, and total angular momentum of 
the ISCO for corotating equal-mass black holes on a 
maximal slice. For the quasiequilibrium data defined in 
this paper, two definitions of the ISCO are listed. One 
uses the minimum in B ADM along sequences where con- 
dition H75fl is satisfied to define the ISCO. The alter- 
native method uses a minimum in Eh along sequences 
where -Mj rr is held fixed. However, we note that the min- 
ima in Eh along sequences that satisfy B75JI are numer- 
ically indistinguishable from the latter. The ISCO for 
the GGB data (listed as HKV-GGB in the table) is de- 
fined as the minima in Eh along a sequence where Mj rr 
remains fixed. Recall that the definition of the ISCO for 
the PN data is consistent with either definition used for 
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the quasiequilibrium data. In addition to the effective 
one-body (EOB) PN data displayed previously in Fig. [HJ 
we also include "standard" PN results for the ISCO as 
reported in Ref. [ilf . 

For corotating quasiequilibrium data, there is very lit- 
tle difference in the results for the two definitions of the 
ISCO. As we will see later, this is not true for the irro- 
tational configurations (see Table ITTTft . In that case, it is 
clear that only the definition in terms of E\> is consistent 
with the PN data. 
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FIG. 10: ISCO configuration for three different choices of 
the lapse boundary condition for equal-mass corotating and 
irrotational black holes computed with the maximal slicin. 
condition. For comparison, results of Refs. [Tftl . I3H Effll 
are included. For post-Newtonian calculations the size of the 
symbol indicates the order, the largest symbol being 3PN. 
"PN standard" pH represents a PN-expansion in the stan- 
dard form without use of the EOB-technique (only 2PN and 
3PN are plotted). 

Figure EH plots binding energy versus orbital angu- 
lar velocity for the ISCO obtained for all three lapse 
boundary conditions for the corotating quasiequilibrium 
data, as well as the corotating results from GGB and 
PN results. All the numerical results are computed on a 
maximal slice. We see that the results for the different 
lapse boundary conditions are essentially indistinguish- 
able. We also see that the PN results converge roughly 
toward the numerical quasiequilibrium results. While 
we would not expect the quasiequilibrium numerical re- 
sults to agree with any of the individual PN results, we 
might expect the GGB result to agree within numeri- 
cal error. All of our numerical results using different 
lapse boundary conditions are essentially indistinguish- 
able. Furthermore, if we use a lapse boundary condition 
that approaches a Dirichlet value of zero on the excision 
surface, the resulting set of boundary conditions is equiv- 
alent to those used by GGB. The difference seen between 



the GGB ISCO and our results may well be due to the 
regularization procedure introduced by GGB. 
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FIG. 11: Constant M\ TT sequence of corotating equal-mass 
black holes. Comparison of post-Newtonian EOB sequences 
with numerical maximal slicing results from HKV and this 
paper. 
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FIG. 12: ISCO configuration for three different choices of 
the lapse boundary condition for equal-mass corotating and 
irrotational black holes computed with the maximal slicing 
condition. Symbols as in Fig. 1101 

There are three rigorously defined gauge-invariant 
global quantities associated with a black-hole binary sys- 
tem: the ADM energy -E ADM , the total angular momen- 
tum J, and the orbital angular velocity as seen at infinity 
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£!o- Figures 151 and ITUl plotted the binding energy E^ (di- 
rectly related to the ADM energy) as a function of J. 
Figures and El plot Eh as a function of fio for the 
same set of sequences and for the ISCO. And, Figs. El 
and El pl°t J as a function of H,q for the same set of 
sequences and for the ISCO. 
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FIG. 13: Constant M- 1TT sequence of corotating equal-mass 
black holes. Comparison of post-Newtonian EOB sequences 
with numerical maximal slicing results from HKV and this 
paper. 



2. Irrotational binary systems 

For the case of irrotational black holes, Fig. 1151 shows 
the binding energy Ef, of the binary system as a function 
of J for all the lapse boundary conditions. As with the 
corotating black holes, we again see that the choice of 
the lapse boundary condition has very little effect on the 
sequence. In the inset to the figure we see that the effect 
is largest at small separations and that the differences due 
to varying the lapse boundary condition are somewhat 
larger than in the corotating case. 
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FIG. 14: ISCO configuration for three different choices of 
the lapse boundary condition for equal-mass corotating and 
irrotational black holes computed with the maximal slicing 
condition. Symbols as in Fig. 1101 



FIG. 15: Constant Mi IT sequence of irrotational equal-mass 
black holes. Maximal slicing is used in these cases, and three 
different excision boundary conditions for the lapse were used. 

Perhaps the most striking difference between the coro- 
tating and irrotational sequences is that the extrema in 
E}, versus J are much less "sharp" in the irrotational se- 
quences. If we consider sequences, either corotating or 
irrotational, that are normalized so that di? ADM = £7odJ 
is satisfied, then the extremum in E ADM will necessarily 
coincide with the extremum in J leading to a very sharp 
cusp in a plot of these quantities. However, because M- lrr 
is not necessarily constant along a sequence with this 
normalization, the extremum in E^ will not necessarily 
coincide with that of J. Thus, we should certainly not 
expect to see a sharp cusp in either Fig. or Fig. ITol 

Another difference between the corotating and irrota- 
tional sequences can be seen if Fig. Here we see that, 
if we demand that the thermodynamic identity (|75|l be 
satisfied along the sequence, then the variation in M- m is 
20 times larger in the irrotational sequences than in the 
corotating sequences. For the corotating sequences, the 
variation in Mj rr due to differences in the lapse bound- 
ary condition was comparable to the average variation. 
For the irrotational sequences, the effect of different lapse 
boundary conditions is clearly negligible. Furthermore, 
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we note that Mj rr is decreasing as the binary separation 
decreases. This behavior is unphysical, as the irreducible 
mass never decreases; therefore, Mi rr should also not de- 
crease during the insipiral of a binary black hole. 
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FIG. 16: Constant Mi rr sequence of irrotational equal-mass 
black holes. Comparison of post-Newtonian EOB sequences 
with numerical maximal slicing results from Conformal Image 
and this paper. 

Figure ITS1 shows a comparison of our irrotational data 
to the effective one-body post-Newtonian results for ir- 
rotational holes as reported in Ref. 35]. First, second, 
and third post-Newtonian results are displayed (labeled 
IR:EOB-lPN, IR:EOB-2PN, and IR:EOB-3PN respec- 
tively in figure legends). The 3PN results correspond 
to the approach labeled "3PN corot. A(u, 0)" in Ta- 
ble I of Ref. [3f|. Also plotted in this figure is the 
first sequence of numerical initial-data solutions for an 
equal-mass black- hole binary in quasicircular orbit, ob- 
tained from inversion-symmetric initial data using an 
effective potential approach|3!j (labeled IR:Conf. Imag- 
ing/Eff. Pot. or IVP-conf in figure legends). 

Again, there is good agreement between all of the re- 
sults at large separation. Also, it appears that the PN 
results are converging toward the irrotational quasiequi- 
librium numerical results, even when the black holes 
are quite close to each other as seen in the figure's in- 
set. We also see that the early numerical results ob- 
tained from the conformal-imaging data 39J differ only 
slightly from the quasiequilibrium results up the location 
of the 1SCO in the quasiequilibrium sequence. However, 
the conformal-imaging sequence extends to much smaller 
separations before encountering its ISCO. 

Table lTTTl displavs the ISCO data for irrotational equal- 
mass black holes on a maximal slice. Again, the ISCO is 
defined in two ways for the QE-sequence data. However, 
unlike the corotating case, the two definitions disagree 



TABLE III: Parameters of the ISCO configuration for irro- 
tational equal-mass black holes computed with the maximal 
slicing condition. Results are given for three different choices 
of the lapse boundary condition and two choices for the defi- 
nition of the location of the ISCO. Layout as in Table ITT1 
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aip = i ADM 


0.148 
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E b 


0.103 


-0.0181 


0.765 


d(a^) _ aV- ADM 

ar lr 


0.145 


-0.0146 


0.761 


E b 


0.101 


-0.0181 


0.766 


Conf. Imag. 


0.166 


-0.0225 


0.744 


1PN EOB 


0.0692 


-0.0144 


0.866 


2PN EOB 


0.0732 


-0.0150 


0.852 


3PN EOB 


0.0882 


-0.0167 


0.820 


1PN standard 


0.5224 


-0.0405 


0.621 


2PN standard 


0.1371 


-0.0199 


0.779 


3PN standard 


0.1287 


-0.0193 


0.786 



dramatically. It is clear that the ISCO, when defined as 
a minimum of Eb along a sequence where M\„ is held 
fixed, is consistent with the results from the PN data. 
This is fortunate given the fact that M- lrr decreases as 
the binary separation decreases along sequences where 
Eq. (|75l) is satisfied. However, it is unclear why these 
E ADM defined ISCOs disagree so significantly from the 
PN data. The Eb-dctincd ISCO data for the various ir- 
rotational sequences is plotted with the corotating data 
in Figs. UUI E3 and IUI 
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FIG. 17: Constant Mi TT sequence of irrotational equal-mass 
black holes. Comparison of post-Newtonian EOB sequences 
with numerical maximal slicing results from Conformal Image 
and this paper. 

Figures ITrjl and UUI plotted the binding energy Eb as a 
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function of J for the irrotational sequences. Figures IT71 
and El plot Ef, as a function of flo for the same set of 
sequences and for the ISCO. And, Figs . 1 1 81 and 1 1 41 plot J 
as a function of f2o for the same set of sequences and for 
the ISCO. 
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FIG. 18: Constant Mi rr sequence of irrotational equal-mass 
black holes. Comparison of post-Newtonian EOB sequences 
with numerical maximal slicing results from Conformal Image 
and this paper. 



A convenient choice for a non-maximal slicing is one 
based on a stationary black hole in Kerr-Schild coordi- 
nates. For the case of a non-charged, non-spinning black 
hole these are also referred to as ingoing Eddington- 
Finkclstein coordinates (cf. Ref. 42] ). The spatial met- 
ric is given by Eq. (|55|l . and after the coordinate trans- 
formation H5b|) it becomes conformally flat. The trace 
of the extrinsic curvature for the Kerr-Schild slicing of 
Schwarzschild is given by Eq. (|58c|) . For a binary system, 
we use a linear combination of two copies of Eq. H58cfl . 
each centered at the location of one of the black holes, to 
define K . 
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B. Non- maximal slicing 

For the case of sequences of corotating or irrotational 
QE initial data obtained on maximal slices (i.e. A" = 0), 
the choice of the lapse boundary condition seemed to 
have very little effect on gauge invariant quantities. This 
lends support to our assertion that the choice of the lapse 
boundary condition is part of the initial temporal gauge 
freedom. To further test this assertion, we should con- 
sider varying other aspects of the freely specifiable data. 
Quasiequilibrium considerations demand that we choose 
Uij = and d t K = 0. This leaves us the options of 
choosing a non-flat conformal metric or a non-maximal 
slice. 

While changing either can affect the content of the dy- 
namical degrees of freedom in the initial data, the confor- 
mal metric is more closely tied to these dynamical degrees 
of freedom and to the spatial gauge freedom. The choice 
of the trace of the extrinsic curvature is usually thought 
of as fixing the initial temporal gauge freedom, although 
as we have seen from the example of the family of maxi- 
mal slices of Schwarzschild in Sec. lIII 01 this is not always 
sufficient to fix this aspect of the gauge freedom. In any 
case, it seems reasonable that the best choice is to vary 
K. 



FIG. 19: Constant Mi rr sequence of corotating equal-mass 
black holes. Eddington-Finkelstein slicing is used in these 
cases, and three different excision boundary conditions for 
the lapse were used. 

Figures IT§1 and I2U1 display the results for both corotat- 
ing and irrotational sequences of equal-mass black-hole 
binaries based on the Kerr-Schild-likc slicing described 
above. In both cases, the same three lapse boundary 
conditions used for the maximal slicing solutions were 
again used. When the holes are at large separation, the 
different lapse boundary conditions cause little variation 
in the results. However, when the holes are close to- 
gether, the different lapse boundary conditions cause sig- 
nificant variation in the sequences. From this example 
it seems that the choice of the lapse boundary condition 
may have a significant effect on QE solutions of the con- 
formal thin-sandwich equations. However, this example 
may be somewhat misleading. 

Maximal slicing (K — 0) is based on a global geomet- 
ric concept that does not depend on the separation of 
the black holes in a binary. For an isolated black hole, 
the Kerr-Schild slicing also has a geometric interpreta- 
tion. However, a linear combination of the traces of the 
extrinsic curvatures for individual black holes does not 
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FIG. 20: Constant Mi rr sequence of irrotational equal-mass 
black holes. Eddington-Finkelstein slicing is used in these 
cases, and three different excision boundary conditions for 
the lapse were used. 



retain this geometrical meaning. Thus, as we vary the 
separation of the black holes in the non-maximal slicing 
sequences, we are also effectively varying the slicing con- 
dition. This effect is weak when the holes are at large 
separation, but becomes significant when the black holes 
are close together. 

In constructing meaningful sequences, everything in 
the construction of the individual models should be held 
fixed except for the separation. For maximal slicing, the 
various choices of the lapse boundary condition choose a 
particular slice from among a family of maximal slices. 
However in the of the Kerr-Schild based slicing, is seems 
likely that the functional form of if as a function of sep- 
aration, and the form of the lapse boundary condition 
conspire to define a different slicing condition for each 
initial-data model considered. Therefore, it seems rea- 
sonable to question the validity of these non-maximal 
slicing sequences. More importantly, we should be cau- 
tious in attributing undue significance to the choice of 
the lapse boundary condition based on this example. 



VI. DISCUSSION 

In this paper, we have refined the QE boundary con- 
ditions defined originally in Ref. and explored both 
single and binary black-hole configurations. The orig- 
inal motivation in deriving these boundary conditions 
was to provide conditions that would be consistent with 
quasiequilibrium configurations. The binary black hole 
initial-data sets constructed in Sec. El are intended to be 
in quasiequilibrium. While the individual black holes in 



a binary cannot be in true equilibrium, it would be useful 
to determine if they are roughly in equilibrium. 

One measure of this is to see how well the QE lapse 
condition (|51|l is satisfied. For stationary black holes, this 
equation holds. However, in the general case, it does not. 
Equation (|51fl defines asc a s the error in this boundary 
condition when applied on the excision boundary of one 
of the holes. 







u 

S B 0.0004 



6/J 
> 
< 



-0.0008 



0.008 



T 




CO: MS - d(ca|r)/dr = 
CO: MS - ay = 1/2 
CO: MS - d(ay)/dr = (ay)/2r 




FIG. 21: Plot of the residual of the quasiequilibrium boundary 
condition on the lapse 1511 of one black hole. Three different 
choices for the lapse boundary condition are shown for coro- 
tating black holes in an equal-mass binary. The upper plot 
shows the average of the residual over the boundary surface. 
The lower plot shows the I/2-norm of the residual. The verti- 
cal dashed line shows the approximate location of the ISCO 
defined by a minimum in Eb- The vertical dotted line shows 
the approximate location of the ISCO defined by a minimum 
in Badm- 

Figure shows both the average value of asc an d 
the Li norm of olbc a s a function of f2o for a corotating 
equal-mass binary. Figure l221 shows the same information 
for an irrotational equal-mass binary. At large separa- 
tions (small fio), we should expect that each black hole 
is nearly in equilibrium. For the corotating sequences, 
mQo rs 0.01 corresponds to a proper separation between 
the horizons of approximately 20m. At this separation, 
\(*bc\l 3 ~ 0.0003. At the ISCO separation, mft ~ 0.11 
corresponding to a proper separation between the hori- 
zons of approximately 4.5m, and |asc|L 2 nas increased 
by a factor of approximately 20. As we might expect, 
the level of violation of the QE lapse boundary condition 
increases steadily as the separation between the holes de- 
creases. However, there is no dramatic increase in the 
violation near the ISCO. For the irrotational sequences, 
\&bc\l 2 begins at mflg w 0.01 at a level approximately 
twice as large as that of the corotating sequence. Near 
the ISCO, it has increased by a factor of approximately 
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FIG. 22: Plot of the residual of the quasiequilibrium boundary 
condition on the lapse j^J of one black hole. Three different 
choices for the lapse boundary condition are shown for irro- 
tational black holes in an equal-mass binary. The upper plot 
shows the average of the residual over the boundary surface. 
The lower plot shows the I/2-norm of the residual. The verti- 
cal dashed line shows the approximate location of the ISCO 
defined by a minimum in Eb- 

It seems that the rate of increase in the violation of 
the QE lapse boundary condition for the irrotational 
sequences is faster than that seen in the corotating se- 
quences. This is not too surprising when we recall that a 
true stationary binary configuration can only be achieved 
for corotating binaries |3r3|. As with the corotating se- 
quences, the level of violation of the QE lapse boundary 
condition increases steadily as the separation between the 
holes decreases, and there is no dramatic increase in the 
violation near the ISCO. 

Another indicator of whether or not each black hole in 
the binary is in equilibrium is given by the value of dtip 
as evaluated on the apparent horizon. We can express 
the time derivative of the conformal factor on any closed 
surface as 



(80) 



\h kl u kl + V20 - (a - /9J H 



Clearly, when the QE boundary conditions in Eqs. 1171) 
and i|39[) are imposed on the excision surface, the last two 
terms in Eq. Ij80|l vanish. Furthermore, in constructing 
QE configurations, we have also demanded that Uij — 
globally. Therefore, the only terms that are possibly non- 
zero on the excision surface are those that involve (3^ . 

For corotating binaries, (3 n l = and we find that the 
time derivative of the conformal factor vanishes identi- 



cally on the excision surface. This is confirmed in our 
numerical results as shown in the upper half of Fig. 1231 
There, we see that dtln?p = to roundoff error. For 
irrotational binaries, the QE conditions require that we 
take /3|| 1 proportional to a conformal Killing vector of hij . 
This implies that /3 y 1 will also be a conformal Killing vec- 
tor of hij. Unfortunately, the operator acting on (3^ in 
Eq. H80f) is not the conformal Killing operator, but rather 



(81) 



While this would vanish if (3 n were a Killing vector of hij , 
it will not vanish if /3y k is only a conformal Killing vector 
of hij , as it will be unless the configuration is truly sta- 
tionary. Again, this is confirmed in our numerical results 
where we find that dthiip ~ 10 -5 when m£l ~ 0.01 and 
grows monotonically as the binary approaches the ISCO. 
These results are shown in the lower half of Fig. 1231 
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FIG. 23: Plot of the 1,2-norm of dt In ip as evaluated on the ex- 
cision boundary of one black hole. The upper half of the figure 
shows the results for corotating black holes in an equal-mass 
binary, while the lower half shows the results for irrotational 
black holes. For both cases, three different choices for the 
lapse boundary condition are shown. In the upper half of the 
figure, the vertical dashed line shows the approximate loca- 
tion of the ISCO defined by a minimum in Eb- The vertical 
dotted line shows the approximate location of the ISCO de- 
fined by a minimum in E A um- In the lower half of the figure, 
the vertical dashed line shows the approximate location of the 
ISCO defined by a minimum in Eb. 

The quasiequilibrium boundary conditions we have de- 
rived and tested in this paper are extremely general. 
Within the conformal thin-sandwich approach, they will 
work for any number of black holes that are to be consid- 
ered in quasiequilibrium, or "isolated" . In this paper, we 
have used several different choices for JC but maintained 
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the assumption of conformal flatness. We emphasize that 
this is not a limitation of the boundary conditions which 
can, in fact, be used with any viable choice for the con- 
formal thrcc-gcomctry specified by 7^ . Furthermore, for 
binary systems, we have only considered the special cases 
of corotating and irrotational black holes. Again, this is 
not a limitation of the boundary conditions which can, 
in principle, produce any desired spin on the individual 
black holes. 

It has been pointed out that the boundary condi- 
tions we have derived are precisely those required to 
construct a black hole satisfying the isolated-horizon 
conditions [T^. I20L I2H Elij . This is not surprising since 
the physical notions underlying an isolated horizon and 
a black hole in quasiequilibrium are essentially the same 
thing. It seems likely that the unified approach offered 
by the isolated-horizon framework will prove useful in 
further understanding the physical content of the binary 
black hole initial data constructed with the quasiequilib- 
rium boundary conditions and to further understand the 
role of the lapse boundary condition. In fact, during the 
final stages of the preparation of this manuscript, we be- 
came aware of a paper by Jaramillo et al.0] that makes 
the connection between our quasiequilibrium boundary 
conditions and isolated horizons more precise. This pa- 
per argues that the lapse boundary condition, Eq. I|51|) 
derived previously in Ref. Q, could be problematic, as we 
have found and discussed, and shows that weakly isolated 
horizon considerations do not restrict the lapse bound- 
ary condition when constructing initial data, consistent 
with our findings here. Furthermore, they suggest an 
alternate boundary condition on the lapse, based on a 
Lie-derivative along the null-generators of the horizon. 
It is not immediately clear that this proposed bound- 
ary condition can work. For the case of a static black 
hole, we have shown that essentially any boundary con- 
dition on the lapse, when combined with the quasiequilib- 
rium boundary conditions, will yield a valid static slice of 
Schwarzschild. It may well be that this proposed bound- 
ary condition is degenerate similar to Eq. I|51[) (cf. our 
discussion in Sec. IIII C|l . 

Clearly, additional work is required to fully understand 
the boundary conditions we have derived, and in particu- 
lar the proper role of the lapse boundary condition. How- 
ever, it is also clear that obtaining appropriate bound- 
ary conditions is not the final issue in the quest to con- 
struct astrophysically realistic binary black hole initial 
data. The most pressing issue is the question of how to 
make a realistic choice for the conformal three-geometry. 
While the errors introduced by the assumption of con- 
formal flatness are not "grave", it is clear that we must 
find a way to allow the physics to dictate the conformal 
three-geometry instead of choosing it a priori. The ap- 
proach al ong these directions outlined in Ref. (see 
also Ref. 46]) is clearly promising. 
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APPENDIX A: COROTATING SEQUENCE 

In this Appendix, we list the numerical results for coro- 
tating equal-mass black holes assuming conformal flat- 
ness, maximal slicing, and using Eq. (|59b|) for the lapse 
boundary condition on both excision surfaces. The data 
has been scaled so that the sequence satisfies Eq. (|75f) 
by following the procedure outlined in Eqs. (|76|) - (|79|l . 
In order to maintain accuracy in the scaling, the max- 
imum coordinate separation between successive models 
was Ad = 0.05. Data in the given tables can be easily 
rescaled to construct sequences with Mj rr held constant. 

In Table ITVl d is the coordinate separation of the cen- 
ters of the excised regions. M m is the irreducible mass 
associated with one of the black holes. E ADM is the ADM 
energy of the system. Qq is the orbital angular velocity of 
the binary system as measured at infinity. is the bind- 
ing energy of the system defined as E\, = £ ADM — 2M- 1IT . 
J ADM is the total ADM angular momentum of the binary 
system as measured at infinity. Finally, £ is the minimum 
proper separation between the two excision surfaces as 
measured on the initial-data slice. 



APPENDIX B: IRROTATIONAL SEQUENCE 

In this Appendix, we list the numerical results for irro- 
tational equal-mass black holes assuming conformal flat- 
ness, maximal slicing, and using Eq. (|59b|) for the lapse 
boundary condition on both excision surfaces. The data 
has been scaled so that the sequence satisfies Eq. (|75f) 
by following the procedure outlined in Eqs. (|76() - (|79|l . 
In order to maintain accuracy in the scaling, the max- 
imum coordinate separation between successive models 
was Ad = 0.05. Data in the given tables can be easily 
rescaled to construct sequences with Mj rr held constant. 

In Table [V] d is the coordinate separation of the cen- 
ters of the excised regions. M- m is the irreducible mass 
associated with one of the black holes. E ADM is the ADM 
energy of the system. SIq is the orbital angular velocity of 
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TABLE IV: Sequence of corotating equal-mass black holes on 
a maximal slice. The length scale is set so that the ADM 
mass of the binary at infinite separation is 1. The ISCO is at 
separation d = 8.28. 



d 



1 



o 



Ja 



40 
35 
30 
25 
20 
19 
18 
17 
16 



15 
14.5 
14 
13.5 
13 
12.5 
12 
11.5 
11 
10.5 
10 

9.5 
9 
8.9 
8.8 
8.7 

8.6 
8.5 
8.4 
8.35 
8.3 
8.28 
8.2 
8.1 
8 



5000000 
5000001 
5000003 
5000006 
5000013 
5000016 
5000020 
5000024 
5000030 
5000039 
5000044 
5000050 
5000059 
5000065 
5000076 
5000088 
5000104 
5000123 
5000147 
5000177 
5000217 
5000268 
5000280 
5000293 
5000306 
5000321 
5000336 
5000352 
5000361 
5000369 
5000373 
5000388 
5000407 
5000428 



-0.0058296 
-0.0065815 
■0.0075478 
-0.0088277 
-0.0105789 
-0.0110046 
-0.0114600 
-0.0119466 
-0.0124654 
-0.0130163 
-0.0133032 
-0.0135969 
-0.0138966 
■0.0142009 
■0.0145079 
■0.0148151 
-0.0151187 
-0.0154139 
-0.0156939 
-0.0159493 
-0.0161675 
-0.0163311 
-0.0163553 
-0.0163761 
-0.0163934 
-0.0164068 
-0.0164160 
-0.0164209 
-0.0164215 
-0.0164209 
-0.0164203 
-0.0164159 
■0.0164054 
-0.0163890 



0.01090 

0.01327 

0.01665 

0.02175 

0.03012 

0.03246 

0.03511 

0.03814 

0.04164 

0.04571 

0.04800 

0.05049 

0.05320 

0.05617 

0.05942 

0.06300 

0.06696 

0.07135 

0.07625 

0.08174 

0.08792 

0.09492 

0.09644 

0.09799 

0.09959 

0.1012 

0.1029 

0.1047 

0.1055 

0.1064 

0.1068 

0.1083 

0.1102 

0.1121 



-0.0058296 
-0.0065816 
■0.0075483 
■0.0088289 
■0.0105816 
-0.0110079 
-0.0114639 
-0.0119514 
-0.0124715 
-0.0130240 
-0.0133119 
-0.0136068 
-0.0139080 
-0.0142140 
-0.0145231 
-0.0148328 
■0.0151395 
■0.0154385 
■0.0157232 
-0.0159847 
-0.0162108 
-0.0163846 
■0.0164112 
-0.0164346 
-0.0164546 
-0.0164709 
-0.0164832 
-0.0164913 
-0.0164936 
-0.0164948 
-0.0164949 
-0.0164934 
-0.0164869 
-0.0164747 



1.2280 
1.1655 
1.1005 
1.0332 
0.9647 
0.9511 
0.9376 
0.9243 
0.9113 
0.8986 
0.8925 
0.8865 
0.8807 
0.8752 
0.8699 
0.8649 
0.8602 
0.8559 
0.8521 
0.8489 
0.8463 
0.8445 
0.8442 
0.8440 
0.8438 
0.8437 
0.8436 
0.8436 
0.8436 
0.8436 
0.8436 
0.8436 
0.8437 
0.8439 



21.81 
19.17 
16.52 
13.83 
11.11 
10.56 
10.00 
9.444 
8.882 
8.316 
8.031 
7.745 
7.458 
7.169 
6.879 
6.587 
6.293 
5.997 
5.699 
5.398 
5.095 
4.788 
4.726 
4.665 
4.603 
4.540 
4.478 
4.416 
4.384 
4.353 
4.341 
4.290 
4.227 
4.164 



the binary system as measured at infinity. is the bind- 
ing energy of the system defined as Ef, = E ADM — 2M m . 
J ADM is the total ADM angular momentum of the binary 
system as measured at infinity. Finally, I is the minimum 
proper separation between the two excision surfaces as 
measured on the initial-data slice. 
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TABLE V: Sequence of irrotational equal-mass black holes 
on a maximal slice. The length scale is set so that the ADM 
mass of the binary at infinite separation is 1. The ISCO is at 
separation d = 8.69. 



d 



E A 



1 



o 



E b 



Ja 



40 
35 
30 
25 
20 
19 
18 
17 
16 
15 



14.5 
14 
13.5 
13 
12.5 
12 
11.5 
11 
10.5 
10 



9.5 
9 

8.9 
8.8 



8.69 
8.6 



8.5 

8 

7.5 
7 

6.9 
6.86 
6.7 
6.5 



5000000 
4999988 
4999963 
4999903 
4999730 
4999661 
4999571 
4999451 
4999286 
4999057 
4998908 
4998729 
4998513 
4998248 
4997922 
4997517 
4997006 
4996356 
4995517 
4994416 
4992944 
4990934 
4990447 
4989924 
4989306 
4988763 
4988117 
4984034 
4977854 
4967902 
4965184 
4964004 
4958660 
4950163 



-0.0058830 
-0.0066618 
■0.0076750 
■0.0090444 
■0.0109891 
-0.0114788 
■0.0120118 
-0.0125938 
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-0.0139306 
-0.0143062 
-0.0147003 
■0.0151140 
■0.0155484 
-0.0160042 
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-0.0186189 
-0.0192006 
-0.0197901 
-0.0199079 
-0.0200252 
-0.0201537 
-0.0202582 
-0.0203733 
-0.0209241 
-0.0213930 
-0.0216779 
-0.0216959 
-0.0216980 
■0.0216708 
-0.0215319 



0.01090 

0.01327 

0.01665 

0.02175 

0.03014 

0.03247 

0.03513 

0.03817 

0.04169 

0.04578 

0.04808 

0.05059 

0.05332 

0.05632 

0.05961 

0.06324 

0.06726 

0.07174 

0.07676 

0.08241 

0.08882 

0.09616 

0.09776 

0.09940 

0.1013 

0.1029 

0.1047 

0.1146 

0.1265 

0.1411 

0.1445 

0.1459 

0.1518 

0.1600 



-0.0058830 
-0.0066594 
■0.0076676 
■0.0090250 
■0.0109351 
-0.0114111 
-0.0119261 
-0.0124839 
-0.0130882 
-0.0137419 
-0.0140878 
-0.0144462 
•0.0148166 
-0.0151980 
-0.0155887 
■0.0159857 
-0.0163846 
-0.0167784 
■0.0171564 
-0.0175020 
■0.0177893 
■0.0179770 
-0.0179972 
-0.0180101 
-0.0180149 
■0.0180108 
-0.0179966 
-0.0177308 
-0.0169638 
-0.0152582 
-0.0147327 
-0.0144988 
-0.0134029 
-0.0115644 



1.2175 
1.1527 
1.0846 
1.0127 
0.9367 
0.9211 
0.9053 
0.8894 
0.8734 
0.8574 
0.8494 
0.8414 
0.8334 
0.8255 
0.8176 
0.8098 
0.8022 
0.7946 
0.7873 
0.7802 
0.7733 
0.7670 
0.7658 
0.7646 
0.7633 
0.7623 
0.7611 
0.7561 
0.7522 
0.7500 
0.7499 
0.7499 
0.7501 
0.7510 



21.81 
19.18 
16.52 
13.84 
11.12 
10.57 
10.01 
9.454 
8.893 
8.329 
8.045 
7.759 
7.473 
7.185 
6.895 
6.604 
6.311 
6.016 
5.718 
5.419 
5.116 
4.811 
4.749 
4.688 
4.620 
4.564 
4.502 
4.189 
3.872 
3.548 
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3.457 
3.351 
3.218 
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